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Super Soft Source (SSS) spectra are powered by nuclear burning on the surface of a white dwarf. The released energy 
causes a radiatively-driven wind that leads to a radially extended atmosphere around the white dwarf. Significant blue 
shifts in photospheric absorption lines are found in the spectra of novae during their SSS phase, being an evidence of 
continued mass loss in this phase. We present spherically symmetric PHOENIX models that account for the expansion 
of the ejecta. A comparison to a plane parallel, hydrostatic atmosphere model demonstrates that the mass loss can have 
a significant impact on the model spectra. The dynamic model yields less pronounced absorption edges, and harder X- 
ray spectra are the result. Therefore, lower effective temperatures are needed to explain the observed spectra. Although 
both types of models are yet to be fine-tuned in order to accurately determine best fit parameters, the implications on the 
chemical abundances are going in opposite directions. With the expanding models the requirement for strong depletion of 
the crucial elements that cause these edges is now avoidable. 
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1 Introduction 

SSS are commonly accepted to be supported by nuclear fu- 
sion on the surface of a white dwarf, where the hydrogen- 
rich burning material is provided by accretion from a com- 
panion star (Kahabka & van den Heuvel 1997). Observed 
SSS spectra show the direct signatures of nuclear burning 
(Bremsstrahlung spectra of ~ 10— 30keV). The white dwarf 
is surrounded by an optically thick atmosphere that emits in 
soft X-rays. In addition to the so-called persistent SSSs like 
Cal 83, some classical novae were observed to go through 
a phase of SSS emission, when their spectrum resembles 
that of Cal 83. This is not surprising, as it is expected that 
the opacity of the ejecta is changing with time, which leads 
to a shift of the peak of the spectral energy distribution 
from the optical via the UV to the soft X-rays (see e.g., 



Gallagher & Starrfieldll 19781) . The X-ray emitting phase is 



reached when the radius of the photosphere has receded 
into sufficiently hot layers. This is a consequence of de- 
clining mass-loss rate which leads to a decrease of the den- 
sity in the envelope. Novae have been observed with ex- 
tremely bright SSS emission, e.g. V1974Cyg by (Krautter 
et al. (1996), much brighter than the persistent SSSs. The 
brightness of nova SSS spectra makes them ideal targets 
for spectral analysis. In this paper we focus on the novae 
RS Oph and V2491Cyg. 

The X-ray gratings on board XMM-Newton and Chan- 
dra yield a spectral resolution that is two orders of mag- 
nitude higher than precursor X-ray satellites ROSAT and 
BeppoSax. The first X-ray grating spectrum of a nova in 



outburst, presented by iNess et alj (12003b . yields deep ab- 
sorption lines. These high quality spectra are oftentimes in- 
terpr eted with the help of h ydrostatic atmosphe re models, 
e.g. dRauch et alj|2005l) and dNelson et al.ll2008l) . But even 
if, in some cases, atmosphere models yield better represen- 
tation of the data than blackbody fits, a good fit alone is not 
a sufficient justification for the model to be realistic. 

An important achievement from the grating spectrom- 
eters is the detection of line shifts in t he SSS spectra of 
novae, as first see n by Ness e t al. (2003) and illustrated by 
Ness et all J20071) and lNessI d2009h . In hydrostatic models 
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these line shifts are not recovered. Blue-shifting the whole 
synthetic spectrum with a fixed amount in order to better fit 
the observation does not physically represent an expanding 
atmosphere. The expansion is a spherical effect in which 
some parts of the atmosphere move towards the observer 
while some other parts move away. Furthermore, the expan- 
sion affects the radial extension of the atmosphere. In this 
paper we present expanding atmosphere models, compare 
them to static models, and show that the results obtained 
from fits to observations are significantly different. 

2 Expanding atmosphere models 

We present first results of new models that were computed 
with the general-purpose, state-of-the-art stellar atmosphere 
code PHOENIX (Hauschildt & Baron 1992, Hauschildt & 
Baron 1999, Hauschildt & Baron 2004). The radiation trans- 
fer equations are solved using the operator splitting method. 
Spherical symmetry is assumed for the atmosphere, and NLTE| 
computation is based on the assumption of statistical equi- 
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librium. PHOEN IX was used to model X-ray novae before 
dPetz et alj|2005b . 



Since 2005, the code has been improved significantly, 
namely, NLTE has been reimplemented with new rates and 
opacities, a new temperature correction method was used, 
broad lines are handled differently, and a new (hybrid, see 
below) atmosphere construct ion method has be en imple- 
mented, based on the idea of lAufdenbergl d2000b . The im- 
proved code is physically more realistic and is a factor 15- 
45 faster. For details we refer to van Rossum (Ph.D. thesis 
in preparation). 

The hybrid atmosphere consists of a hydrostatic base 
with an expanding envelope on top. The hydrostatics are 
parametrised in the classic way by \og(g), ro and T c s . Here, 
ro is the outer radius of the static core. For the expanding 
envelope, a velocity field is defined, and from the conti- 
nuity equation for a constant mass loss rate, M, a density 
profile follows. The velocity fields assumed in the models 
presented here are described by the beta-law (Lamers & 
Cassinelli 1999). 



v(r > ro) = Voo(l - r /r) 



(1) 



The additional parameters that describe the expanding en- 
velope are the terminal velocity and the exponent j3. 

In this paper we present the first results of the models, 
based on which more refined models can be calculated. As 
a start, we assumed solar abundances of the elements H, He, 
C, N, O, Ne, Al, Mg, Si, S, Ar, Ca, and Fe (yielding roughly 
7500 independent levels with 150,000 transitions, all calcu- 
lated in NLTE). We emphasise that the models presented in 
this paper have not been fine tuned, i.e., the full parameter 
space has not been explored yet. 

3 Results 

In Fig. [U we illustrate the effect of the mass loss parame- 
ter on the model spectrum. By increasing the mass-loss rate 
the ionisation edges become weaker and thus the spectra 
become harder. As a consequence, dynamic models that are 
fitted to observed spectra typically yield a lower effective 
temperature than static models. 

In Fig. |2] we show the observation of RS Oph that was 
taken on day 39.7 after outburst (for an observation log we 
refer to Ness 2009 and references therein). A small grid of 
models was calculated with a terminal velocity of v^, = 
1200 k m/s, according to the blue-shifts observed for this 
object dNess et al.l 12007). We achieve the best agreement 
with X" c ff = 5.5 • 10 5 K, logg = 8.4, a wind mass loss 
rate of M = 1.5 • 10~ 8 M /yr, R cS = 7900km, and iV H = 
4.8 • 10 21 cm- 2 . 

Our interstellar extinction model is based on absorption 
cross sections in the X-ray domain compiled by Balucinska- 
Church & McCammon (1992) for 17 astrophysically impor- 
tant species. Specifically, the effective extinction curve was 
calculated with Balucinska-Church & McCammon's code, 
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Fig. 1 Example of the impact of the mass-loss rate on 
the model spectrum. A range of models with T c g = 6.0 • 
10 5 K and — 2400 km/s and different mass-loss rates is 
shown. With increasing mass-loss rate the ionisation edges 
at 25.3 A, 18.6 A and 16.8 A (from the elements C, N and O, 
respectively) become weaker, and the flux between 14 and 
30A increases. 



assuming solar abundances. The total extinction is then pro- 
portional to the hydrogen column density. We would like to 
point out that, owing to the large interstellar column den- 
sity, the soft part of the spectrum is significantly underrep- 
resented when fitting a model to an absorbed spectrum. In 
the top panel of Fig.|2j it can be seen that the emission that 
is produced by the source (thus before traversing through 
the interstellar medium), peaks at much longer wavelengths 
than the observed spectrum (bottom panel in grey shadings; 
see also Fig. 1 in Ness 2009). Since the ultimate goal is 
to reproduce the original (unabsorbed) spectrum (note that 
the interstellar extinction has no physical relation to the ob- 
ject we try to model), reproduction of the soft tail is at least 
as important as the rest of the spectrum. Since the soft tail 
is primarily determined by the interstellar extinction, this 
leads to the conclusion that iVn cannot be treated as a free 
parameter in the fitting procedure. Instead, A^h must be de- 
termined to reproduce the soft tail. This restriction in free 
parameter space grievely limits the possibilities in finding a 
good fit. Accounting for this effect, especially, hydrostatic 
models turn out to poorly fit the observations (see for exam- 
ple the bottom panel in Fig. [3}. The usual goodness criteria 
(e.g., x 2 ), determined from fitting a model to the absorbed 
spectrum only yield a description of how well the data, es- 
pecially the brightest parts of the spectrum, are reproduced. 
But that does not necessarily indicate whether the physi- 
cal description of the object is realistic or not. We therefore 
fit the parameter Nn by eye, considering both the absorbed 
and unabsorbed fit quality. 

In Fig.[5J the results for the more recent nova V249 1 Cyg 
(2008) are presented, where we compare with an XMM- 
Newton observation taken on day 49.7 after outburst (grey 
shadings). In the top panel, the best model that has so far 
been found from a small grid of models with Voo = 4800 
km/s is shown, yielding the parameters T c s = 6.0 • 10 5 K, 
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RS Oph (day 39.7, LETGS) MH=4.8e21 cm" 2 




Wavelength [A] 

Fig. 2 Comparison of a first cut PHOENIX model spectrum with the observation of RS Oph. The binwidth is 0.025 A. 
In the top panel, the observed spectrum is corrected for interstellar absorption, yielding the non-absorbed spectrum. The 
interstellar transmission coefficient is shown as dashed line. In the bottom panel, the model was corrected for absorption. 
The continuum slope is reproduced by the model longwards of 14 A. At shorter wavelengths the flux is dominated by 
sources other than the SSS atmosphere, e.g. emission from the shock. These contributions are observed already before the 
SSS phase starts and are not included in the model. 



log g = 8.6, a wind mass loss rate of M — 2 ■ 10 8 M Q /yr, 
R cS = 5300km, and N H = 2.8 • 10 21 cm" 2 . Good agree- 
ment between observation and model is found, and it can be 
expected to yet improve significantly when non-solar abun- 
dances are used and the atmospheric structure is fine-tuned 
to the observation. In the actual comparison in Fig. [3] many 
lines do not yet accurately match their observed strengths 
and the continuum is too faint for shorter wavelengths. The 
latter point may, apart from the model that lacks tuning, be 
attributed to non-atmospheric X-ray sources that are not in- 
cluded in the model. In the bottom panel, the same obser- 
vation is shown in comparison to a plane parallel and static 
PHOENIX model with no mass loss that has been tuned to a 
similar degree as the expanding model. The temperature in 
the static model is significantly higher, which is necessary 
to compensate the lack of hard emission shortwards of the 
N VII and Ovil absorption edges at 18.6 A and 16.8 A, re- 
spectively. As a consequence, a smaller effective radius is 
needed. Note that the obtained bolometric luminosities are 
the same, which is primarily a result of the constant Nh ob- 
tained for both fits. 

But more interesting are the different implications for 
abundances. For example, in the hydrostatic model the N VII 



edge at 18.6 A is too strong, and also the N VII line at 24.8 A 
and the N VI line at 28.8 A, leading to the interpretation that 
the N abundance must be subsolar. In contrast, in the ex- 
panding model the Nvil edge at 18.6 A is very weak, the 
Nvi edge at 22.5 A is too weak and also the Nvil line at 
24.8 A and the N VI line at 28.8 A, so that according to the 
new type of models the N abundance should become super- 
solar. These discrepancies must be examined in detail with 
the help of fine-tuned models, which are not yet available at 
the time of writing. 

4 Discussion and Conclusions 

The comparison between the hydrostatic and the new dy- 
namic model in Fig. [3] demonstrates the essence of treat- 
ing the mass-loss in the model atmosphere. Obviously, the 
model spectra are sensitive to the details of the atmospheric 
structure. Comparing observations to models that account 
for the mass-loss can lead to very different conclusions than 
comparison with static atmospheres. Since the continued 
mass-loss is undeniably confirmed by the observations, this 
calls for caution with the interpretation of nova X-ray spec- 
tra using static atmospheres. For example, the conclusion 
from fine tuning of a static atmosphere model to observa- 
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V2491 Cyg (day 49.7, RGS) MH = 2.8e21 cm" 




, 7^ = 600 kK 
log 3=8.6 

dM/dt=2.0e-8 A/sun/yr 
•u_inf = 4800 km/s 
, fl=5300 km 
log Lbol = 37.4 erg/s 



25 

Wavelength [A] 

Fig. 3 XMM-Newton RGS spectrum of V2491 Cyg (2008). The binwidth is 0.02 A. Two different models are shown 
(red curves), a radially extended model with mass loss (top panel) and a white dwarf (compact, hydrostatic) model (bottom 
panel). No good fit could be found with white dwarf models. The best temperature is much higher and the radius much 
smaller than for the model with mass loss. One important difference is the absence of strong ionisation edges in the 
dynamical model, like at 25.3 A, 18.6 A and 16.8 A. Note that the models use identical (solar) abundances but give very 
different implications on being either under- or overabundant (see the discussion in the text). 



tions of RS Oph, presented bv lNelson et al.1 ( 120081) . that the 
C/N abundance ratio is 1000 times subsolar can expected to 
be very different when dynamical models are used. 

The models presented in this paper all use solar abun- 
dances. In the SSS phase the observed radiation originates 
from much deeper regions than in earlier stages after the 
outburst. The material in those regions is a mixture of dred- 
ged-up material from the white dwarf, accreted material from| 
the companion, and reaccreted ejecta of the initial outburst. 
Therefore, the chemical composition that is "seen" at this 
stage could be significantly different from other stages. Pre- 
sently, there are no detailed theoretical and observational 
constraints on this mixture. Although solar abundances mi- 
ght be unrealistic, it is still a good standard set to start with. 
Detailed abundance analysis with both static and dynamic 
models will be needed in order to draw definitive conclu- 
sions. This must be performed with simultaneous fine-tu- 
ning of the atmospheric structure, a task that is scheduled 
for the near future. 

The preliminary models presented in this paper are a 
promising start. More sophisticated and better fine-tuned 
models will increase the accuracy of the parameters and the 
reproduction of the observed spectra. 
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